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Abstract 

We examine the charmless hadronic B decay modes in the context of i?-parity 
violating [I^p ) supersymmetry. We try to explain the large branching ratio 
(compared to the Standard Model (SM) prediction) of the decay — > rj'K^. 
There exist data for other observed rj^'^ modes and among these modes, the 
decay B^ rjK*^ is also found to be large compared to the SM prediction. 
We investigate all these modes and find that only two pairs of I^p coupling can 
satisfy the requirements without affecting the other B PP and B VP 
decay modes barring the decay B (j)K. From this analysis, we determine the 
preferred values of the I^p couplings and the effective number of color N^- We 
also calculate the CP asymmetry for the observed decay modes affected by these 
new couplings. 
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I. INTRODUCTION 



For the last few years, different exeperimental groups have been accumulating plenty of data 
for the charmless hadronic B decay modes. The CLEO the Belle and the BaBar 

collaboration [0-0] are providing us with the information on the branching ratio (BR) and 
the CP asymmetry for different decay modes. A clear picture is about to emerge from these 
information. 

Among the B PP {P denotes a pseudoscalar meson) decay modes, the branching ratio 
for the decay rj'K^ is found to be still larger than that expected within the Standard 

Model (SM). The SM contribution is about 3a smaller than the experimental world average 
(see Fig.l). Among the B VP {V denotes a vector meson) decay modes, the experimentally 
observed BR for the decay B^ rjK*^ has been aloso found to be 2cr larger than the SM. 
The decay B^ (pK^ has been observed recently, and the BR for the newly observed decay 
B"^ — > rjK*^ is also now available. 

In this paper, we address these large BR problems of B^^^^ — > 7^(')_ft'='=(*o) systems using im- 
parity violating (^p) supersymmetric theories (SUSY). The effects of Ijip couplings on B decays 
have been investigated previously in the literatures |jT3|,n|, where attempts were made to fit 
just the large BR for B^ rj'K"^ |T^. At present, we have many more available results. Some 
of these results are concerned with decay modes involving r/*^'^ and these modes are influenced 
by the same coupling that affects B^ —>■ rj'K^. For example, the decay modes -B^ —>■ rjK*^, 
B^ — i> 7]K*^, B^ rj'K^ are affected by the new couplings which cure the large BR problem 
of B"^ rj'K"^. Hence, it is natural to investigate these newly observed decay modes and 
try to see whether all the available data can be explained simultaneously. We also need to be 
concerned about the other observed (not involving r]'^'^) B PP and B — > VP decay modes, 
which could be influenced by these new couplings. Our effort is not to affect the other modes 
as much as possible, since except for B r]^'^K^*'> decay modes, the other observed modes fit 
the available data well [^3|,|TB| within the SM. Further, using the preferred values of different 
parameters (e.g., new couplings etc.), we also make predictions for CP asymmetrey for these 
observed modes which will be verified in the near future. 

We organize this letter as follows. In section II, we give a very brief introduction to the SM 
and Hamiltonian, and list the possible operators that can contribute to charmless B 
decays. We discuss the B PP and B PV decay modes in section III. The new physics 
contributions to different decay modes are also discussed. In section IV, we show how ^ 
can explain the branching ratio of B ^ t]^''^K^*'^ decay modes without jeopardizing many other 
B PP and B VP decay modes. We also discuss the CP asymmetry of these decay modes. 
We conclude in section V. 
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II. EFFECTIVE HAMILTONIAN FOR CHARMLESS DECAYS 



The effective Hamiltonian for cliarmless nonleptonic B decays can be written as 

IVubV:^ E - VttV,l E + h.c. (1) 



V2 



i=l,2 i=3 



The Wilson coefficients (WCs), q, contain the short- distance QCD corrections. We find all 
our expressions in terms of the effective WCs and refer the reader to the papers |17-^ for a 



detailed discussion. We use the effective WCs for the processes b sqq' and h — >■ dqq' from Ref. 



| T9|| . The regularization scale is taken to be = mfe. In our discussion, we will neglect small 
effects of the electromagnetic moment operator O12, but will take into account effects from the 
four-fermion operators Oi — Oio as well as the chromomagnetic operator On. 

The Iflp part of the superpotential of the minimal supersymmetric standard model (MSSM) 
can contain terms of the form 

= K,UH2 + K,tL^L,El + K,,L,Q,Dl + K,kU-DPl (2) 

where Ei, Ui and Di are respectively the i-th type of lepton, up-quark and down-quark singlet 
superfields, Li and Qi are the SU(2)i doublet lepton and quark superfields, and H2 is the Higgs 
doublet with the appropriate hypercharge. From the symmetry reason, we need A^-fc = —Xjik 
and X'-ji^ = —X'ikj- The bilinear terms can be rotated away with redefinition of lepton and Higgs 



superfields, but the effect reappears as As, A's and lepton-number violating soft terms pi 
The first three terms of Eq.(|]) violate the lepton number, whereas the fourth term violates the 
baryon number. We do not want all these terms to be present simultaneously due to catastrophic 
rates for proton decay. In order to prevent proton decay, one set needs to be forbidden. 

For our purpose, we will assume only A'— type couplings to be present. Then, the effective 
Hamiltonian for charmless nonleptonic B decay can be written as 



HeffiP djdkdn) = df^.^[dnalLdjl3 dkplfiRba] + d^j^J^dnalL^P dk/Sl^iRdjc 



with 



E 



u 



R 

'jkn 



'^ijk^in3 

3 \/ \'* 

E^ijn^ikS 

i=l °"''eiL 



d^ 



jkn 



E4^, U,k,n = l,2) 



-1 ^^II 



(4) 



{j,k = l, n = 2). 
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where a and j3 are color indices and 7^^ = 7^(1 ± 75). The leading order QCD correction to 
this operator is given by a scaling factor / — 2 for = 200 GeV. 

The available data on low energy processes can be used to impose rather strict constraints on 
many of these couplings |22-Q. Most such bounds have been calculated under the assumption 
of there being only one non-zero coupling. There is no strong argument to have only one 
-^p coupling being nonzero. In fact, a hierarchy of couplings may be naturally obtained ||23| on 
account of the mixings in either of the quark and squark sectors. In this paper, we try to find 
out the values of such couplings for which all available data are simultaneously satisfied. An 
important role will be played by the AJjgj -type couplings, the constraints on which are relatively 
weak. 



III. B ^ PP AND B ^VP DECAY MODES 



We consider next the matrix elements of the various vector (V^) and axial vector {A^) quark 
currents between generic meson states. For the decay constants of a pseudoscalar or a vector 
meson defined through 

(0|^jP(p)) = ifpp^ 
{Q\V^\V{e,p)) = fvmve, , 

we use the followings (all values in MeV): 



(5) 



= 215, fK* = 225, /„ = 215, U = 132, /k = 162, L, = 146, /, 



180. 



(6) 



The decay constants of the mass eigenstates r] and rj' are related to those for the weak eigenstates 
through the relations 



fs 



fl 



/,^ = ^cos^--|sin^. 



2-^ cosO — ^ sin 6*. 
V6 V3 



The mixing angle can be inferred from the data on the 77 decay modes |]25| to be ^ —22°. 
The B ^ P matrix element can be parameterized as 



{P{p')\V,\B{p))=[{p' 

and the B V transition is given by 

iy{e,p')\{y,-A,)\B{p)) 
2V 



— mj 



1 1 

B^P , ~ 



+ 



B^P 



(7) 



niB + mv 



{niB + mv)Ai e* +e* -q I -A2 



rriB + mv 



(8) 



2m 



V- 



{Aq - ^3 
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with 

2mvA3 = {niB + my)Ai - (m^ - mv)A2. 

The quantities Fq{*^ , V and Aq, 1,2,3 are the hadronic form factors and their values are given in 
the next section. 

The $,p part of the amphtude of —>■ rj'K^ decay is 



= (d 

+ {d 



R 
121 



L 

121 



"112 



d^ 
"112 



222 "222 



m 



[rris + mu){'mb - niu) 



(9) 



where ^ = 1/Nc [N^ denotes the effective number of color), m = m^,/(mfe — m^) and 

Am, = {M2\J^\B) (Mi|J,jO), 
Am, = {M,\j'nB) (M2|J/^|0). 

J and J' stand for quark currents and the subscripts h and / indicate whether the current 
involves a h quark or only the light quarks. Analogous expressions hold for — > rjK"^ where 
we have to replace A^, by AJ^, A^, by and m^i by m^. Replacing a pseudoscalar meson by 
a vector meson, we also get similar expressions for the amplitudes of B^^'^'^ 'q'K*'^^^^ modes. 
The part of the amplitude of i? ^ (pK decay mode involves only ^222 and (i^225 ^'^d 



'^999 + C^999 ] [^A^] 



(10) 



where A^ = {K\J^\B) {(/)\Ji^\0). 



IV. RESULTS 



In our calculation, we use the following input for the Cabibbo-Kobayashi-Maskawa (CKM) 
angles: 

iKfel = \Vts\ = 0.04, iKfcl = 0.087|Kfe|, IK.I = \Vcd\ = 0.218, (11) 
\Vud\ = 0.9722, \Vcs\ = 0.9740, \Vtb\ = 0.9988, 
/? = 26°, 7 = 110° (80°). 

We first try to explain the large branching ratio of B"^ rj'K^. The observed BR for this 
mode in three different experiments are 0,0,|1^ 



B{B^ r]'K^) X 10^ = 80+^° ± 7 (CLEO), 70 ± 8 ± 5 (BaBar), 79^}^ ± 9 (Belle). (12) 

The three results are close and we use the world average of them: (75 ±7) x 10~^. The maximum 
BR in SM that we find is 42 x 10"^ (Fig. 1). In the SUSY framework, we find that the 
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positive values of (if22 and negative values of ^222 can increase the BR, keeping most of the other 
B PP and B VP modes unaffected. The other combinations are either not enough 
to increase in the BR or affect too many other modes. We divide our results into two cases; 
Case 1: we use only (i^22 (positive values) and 

Case 2: we use a combination of ^^22 (positive values) and 6/222 (negative values). 

Let us start with Case 1. We first discuss the case of 7 = 110°. In this scenario we use 
rris (at scale) = 85 MeV. In Fig. 1, we plot the BR for the decay B^ t]'K^ as a function 
of ^. We have used IA'3231 = IA'3221 = 0.04, 0.06, 0.08 and m^usy = 200 GeV. We take ^^22 to be 
positive. The large branching ratio can be explained for A' > 0.05. In our calculation, we use 
the following form factors: 

A^^^ = A^^P = = 0.45, ^ = Fo^f''^ = = 0.29, F^ff = 0.26. (13) 
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FIG. 1. The BR for the decay rj'K^ vs ^. The solid line is for the SM. The dashed, dotted 

and dot-dashed lines correspond to IA323I = IA322I = 0.04, 0.06, 0.08, respectively. The bold solid lines 
indicate the experimental world average bound. 
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FIG. 2. The BR for the deacy cpK^ vs ^. The solid line is for the SM. The dashed, dotted 

and dot-dashed lines correspond to IA323I = IA322I = 0.04, 0.06, 0.08, respectively. The bold solid lines 
indicate the experimental world average bound. 
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FIG. 3. The BR for the deacy r]K*^ vs The sohd line is for the SM. The dashed, dotted 



and dot-dashed hnes correspond to IA323I 
indicate the experimental bound. 
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FIG. 4. The BR for the deacy S° rjK*^ vs ^. The solid line is for the SM. The dashed, dotted 



and dot-dashed lines correspond to IA323I 
indicate the experimental bound. 
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0.04, 0.06, 0.08, respectively. The bold sohd hnes 



In Fig. 2, we plot the BR for the deacy (pK^^ for the same set of couplings. The 

observed BR of this mode by the CLEO [0 is (in lO"'') S.St^i ± 0.6. The Belle and the 
BaBar collaboration have also observed this mode RITI]] with BRs (in 10~^) 11.2^2 i 1-4 and 



7.7^|;4 ± 0.8, respectively. From the figure we see that the BR is increasing with |A'|. The 



BR for I A'. 



3231 



I A' 



3221 



0.06 and ^ ~ 0.25 — 0.3 is consistent with the world average bound 



(7 — 9) X 10 ^. Now combining Fig. 1 and Fig. 2, we find that |Ag 



3231 



lAl 



3221 



0.06 is allowed 



by both decay modes for ^ ~ 0.25. We note that the BaBar's number for this mode is quite 
close to the value observed by the CLEO, but the Belle's number is more than 2a away from 
the CLEO's central value. We hope that this discrepancy will be sorted out in the near future. 

In Fig. 3, we exhibit the BR for the deacy — > rjK*^ as a function of ^. The observed BR 
of this mode B is B{B^ r]K*^) x 10^^ = 26.4t|^ ± 3.3. We find that the solution we have 



got from our previous two decay modes holds in this case: i.e., IA323I 
for ^-0.15-0.25. 



I a; 



3221 



0.06 is allowed 
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TABLE I. The branching ratios (B) for B decays into PP and VP final states at ^ = 0.25. 



Decay modes 


B (10-6) 

^ _ 1 1 nO 
y — iiu 


B (10-6) 

— RCp 
J — ou 


Experimental i3 (10 ^) 


B° TT+TT- 


5.96 


6.0 


4.3lli±0.5 (CLEO i), Btlit'^oi (Belle @) 

/I 1 _|_ 1 -1-0 7 /'RpRq-t IIqII'^ 


B+ ^ K+7r° 


9.59 


9.19 


1 n 8+2-1 _|_ 1 n floU 
iu.o_]^ g m i.u 


B+ ^ K^7r+ 


13.6 


15.63 


18.2lt[;±1.6 g, 13.7l™| @ 
18 2+o'i^ ± 2 [01 




15.10 


14.50 


17.2+^-5 ±1.2 m, 19.3+i^-^+ii [1 

' — 2.4 —3.2—0.6 lO 

16.7 ± 1.6 lb 1.3 Hi 


BO ^ K^TT^ 


8.68 


9.87 


14.6+t-?+2-t Pi, 16.0+^-2+2.5 ^ 

"—5.1—3.3 O' —5.9—2.7 O 

8.2i|i ± 1.2 i 


B+ UJTT+ 


9.518 


9.0 


11.3l|^ ± 1.5 (CLEO §), 6.6tli ± 0.7(BaBar 0) 


B+ ^ pO^+ 


9.418 


9.45 


10.4l|^ ± 2.1 §, 24 ± 8 ± 3 


B^ — > p^TT^ (sum) 


35.23 


35.03 


27.6l^:^ ± 4.2 [| 


B+ ^ujh 

BO K* + TT- 


10.87 
5.93 


12.18 
4.15 


14.3l|6+2.i (CLEO i) 



Since all the parameters are fixed, now it is interesting to see whether the decay rjK*^ 
fits in the allowed region. The observed BR for this mode by CLEO collaboration 0] is (in 
10-6) 13. 814^6 ^ The Belle and the BaBar collaboration have also observed this mode P,p!0[] 
with BRs (in IQ-^) 19.8t|^ ± 1.7 and 21.2t|^ ±2.0, respectively. The SM BR is very small and 
cannot explain the experimental data. In Fig. 4, we plot rjK*^, where the world average 

value of the data is expressed as the bold solid fine. We find that the dotted line (|A'| = 0.06) is 
allowed for ^ < 0.2. But, the estimated BR for IA323I = IA322I = 0.06 at ^ ~ 0.25 is just below 
and very close to the lower bound of the average data. In fact, IA323I = IA322I = 0.06 is allowed 
by both the CLEO and the Belle data for ^ ^ 0.25. Only the BaBar's number is a little bit 
larger than our estimated BR at 0.2 < ,^ ^ 0.25. 

There exists a result for another decay mode involving rj', i.e., rj'K^ . For |A'| = 0.06 

and i = 0.25 we find the BR is 88.3 x lO'^. The CLEO bound is (in lO'^) 89t\l ± 9. The 
Belle and the BaBar collaboration also have reported observation of this mode but with less 
significance, and their results are (in 10-6) 55t\l±8 and42tli±4, respectively. We see that the 
results for this mode from the three experiments are not consistent. Our result in this scenario 
can be taken as a prediction which is in agreement with the CLEO result. 
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TABLE II. The branching ratios (B) and the CP asymmetries for B 7]^'^K^*') and B — > <j)K. 





6 = 0, 




6 = Ifr 




6 = 0, 




5 = 55*^ 






7 = 110'^, 




7 = 110^ 




7 = 80^, 




7 = 80*^ 




mode 


B X 10^ 


Acp 


B X 10^ 


Acp 


B X 10^ 


* 

Acp 


B X 10^ 


Acp 


B+ n'K+ 


68.9 


0.01 


68.3 


0.04 


82.1 


0.01 


68.3 


0.11 


B+ r]K*+ 


36.4 


0.03 


36.4 


0.04 


36.5 


0.03 


32.7 


0.09 




88.3 


0.00 


86.8 


0.03 


110.2 


0.00 


87.1 


0.12 


BO _^ ^K*o 


14.0 


-0.39 


14.6 


-0.42 


14.8 


-0.28 


20.4 


-0.56 


B+ (pK+ 


7.11 


0.00 


6.97 


0.04 


7.10 


0.00 


5.76 


0.14 



The B{B — > Xsw) can put bound on A322A3*23 in certain limits. Using the Refs. p6| and the 



experimental limit (BR< 6.4 x 10~^) on the above BR we find that A' < 0.07. However, if 
we go to any realistic scenario, for example grand unified models (with R parity violation), we 
find a natural hierarchy among the sneutrino and squark masses. The squark masses are much 
heavier than the sneutrino masses and the bound does not apply any more. 

The other observed B PP and B — > VP decay modes are listed in Table I for C, = 0.25 
and we find that the BRs are within the experimental limits. Our result on K*^7t~ is 

compatible within 2a range, but this measurement still involves large error. 

We can use smaller value of 7, e.g., 7 = 80'^ to fit the B r]^'^K^*^ data. In this scenario 
we use rUs (at scale) = 75 MeV. In Table II we show the BRs for B rf^'^K^*^ and in Table 
I we show the BRs for the other observed B —>■ PP and B —>■ VP decay modes in the case of 
7 = 80°. Again, we find the fit is reasonable. In this case, we use Aq^^* = 0.4 and keep the 
other inputs unchanged. 

We now calculate the CP asymmetry for different B 'q^'^K'^*^ modes. The CP asymmetry, 
Acp, is defined by 

_ B{B ^ f) - B(B^ f) 

B{B f) + B{B^f)' ^ ' 

where B and / denote a B meson and a generic final state, respectively. Let us define 

-^323-^322 ~ l-^323-^322k ' (15) 

where 5 denotes the phase difference between A323 and A322- So far we have discussed the 5 = 0° 
situation. In Table II, we calculate the CP asymmetries for B rjK modes for different values 
of 6 and 7. The maximum values of 6 allowed by the BR of B^ — > rj'K^ are S = 15° for 
7 = 110° and 6 = 55° for 7 = 80°. We find that Acp is very large for i?° rjK*^ mode and is 
predicted to be -39% (-42%) for 5 = 0° (15°), 7 = 110°, and -28% (-56%) for 5 = 0° (55°), 
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7 = 80°. Acp for B"^ — > (pK^^ is large (14%) for S = 55^ and 7 = 80°. The other modes are 
also found to be large (~ 10%) for the above set of parameters. 

Case 2: We use the same form factors as in the case 7 = 80° of Case 1 and we use 
rus (at ruh scale) = 85 MeV and 7 = 110°. We now use the combination of (i^2 '^222- 
We assume 

^'2i'2i'2 — ^'2'2i2 ' this scenario, the coupling part of the amplitude in 5 — (pK 
decay mode canceled exactly (Eq. 10). (In fact, our solution still works when the cancellation 
is incomplete by about 5%.) But we still have contributions to B"^ —>■ rj'K^ (Eq. 9) and to 
increase the BR we choose (i^22 be positive. There is no ^ contribution to the other B —>■ PP 
and B VP modes in this case as well. 
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FIG. 5. The BR for the deacy rj'K^ vs ^. The solid line is for the SM. The dashed, dotted 

and dot-dashed lines correspond to |A'| = 0.035, 0.052, 0.07, respectively. The bold solid lines indicate 
the experimental bound. 
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0.2 0.4 0.6 0.8 1 

FIG. 6. The BR for the deacy rjK*^ vs ^. The solid line is for the SM. The dashed, dotted 

and dot-dashed lines correspond to |A'| = 0.035, 0.052, 0.07, respectively. The bold solid lines indicate 
the experimental bound. 



In Fig. 5, we plot the BR for the deacy 5^ rj'K^ as a function of ^. We have used 
|A'| = 0.035, 0.052, 0.07, and m^nsj = 200 GeV. In this case, the large branching ratio can be 
explained for A' > 0.05. 
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FIG. 7. The BR for the deacy S° ^ r]K*^ vs ^. The soHd Une is for the SM. The dashed, dotted 
and dot-dashed Unes correspond to |A'| = 0.035, 0.052, 0.07, respectively. The bold solid lines indicate 
the experimental bound. 



TABLE III. The branching ratios {B) and the CP asymmetries for B rj^'^K^*^ and B (pK. 





(5 = 




5 = 20° 




mode 


B X 10^ 


Acp 


B X 10^ 


Acp 


B+ r]'K+ 


69.3 


0.01 


68.0 


0.05 


B+ r]K*+ 


27.9 


0.04 


27.8 


0.05 


B^ ^ ri'K^ 


107.4 


0.00 


104.5 


0.05 


BO _^ ^K*o 


20.5 


-0.71 


21.1 


-0.72 


B+ (l)K+ 


6.56 


0.00 


6.56 


0.00 



In Fig. 6 and Fig. 7, we plot the BRs for —>■ rjK*^ and 5° rjK*^. Combining Figs. 
5, 6 and 7, we find that | A' | =0.052 and ^ ~ 0.4 — 0.6 can explain all the data. In this scenario 
our result on rjK*^ is allowed by the world average bound. The BR of 5° —>■ rj'K^ is 

107.4 X 10-6 for |A'| = 0.052 and ^ = 0.45 and is allowed by the CLEO data. 

The BR for the deacy —>■ (pK^ does not have a contribution due to the cancellation. 
The SM line (solid) in Fig. 2 needs to be used in this case and we find that ^ ^ 0.45 is allowed. 
The BRs of the other observed B PP and B — VP modes do not get affected by the new 
couplings and these modes seem to fit the data reasonably well for ^ ~ 0.3 — 0.5 [|16|. We also 
calculate the Acp for this case. Since we have assumed that (i^22 ~ "'^222 > phase difference 
between A322 and A332 is (5 + vr), where S is the phase difference between A323 and A322- 

In Table III, we calculate the BRs and the CP asymmetries for B 7]^'^K'^*^ and B (pK 
for different values of 6. The maximum value of 6 allowed by the BR for B^ r]'K^ is 5 = 20° 
for 7 = 110°. As in Casel, we find that Acp is very large for 5° —>■ r]K*^ mode and is predicted 
to be -71% (-72%) for 5 = 0° (20°). The other modes are found to be ~ 5% for the above 
set of parameters. 
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V. CONCLUSION 



We have studied B rj^'^K^*^ modes in the context of supersymmetric theories. We 
have isolated the necessary ^ couphngs (A323A322 and A322A332) that satisfy the experimental 
results for the BRs of these modes reasonably well. The Standard Model contribution is less 
than 2(7 for some of these modes. These new couplings do not affect any other B — > PP and 
B — > VP modes except for the decay B^ — > (l)K^. We have shown that the calculated BR for 
B^ — > (t)K^ agrees with the experimental data. 

We found solutions for both large and small values of 7 = llO*^, 80° and two different values of 
1/Nc) ~ 0.25, 0.45 for two different scenarios. For our solutions, we need |A'| ~ 0.05 — 0.06 
for msusy = 200 GeV. Using these preferred values of the parameters, we calculated the CP 
asymmetry of different observed modes affected by the new couplings and found that the CP 
asymmetry of B° — > rjK*^ is large (~ 28% — 72%) and the CP asymmetry of other modes also 
can be around 10 — 15%. 

These new coupling can be also examined in the RUN II through the associated production 
of the lightest chargino (xf ) ^^^^ the second lightest neutralino (x^)- xf X2 decay into lep- 
ton plus 2 jets (for example), where one of the jets is a 6 jet. The final state of this production 
process contains 2 leptons plus 4 jets. So the signal is quite interesting and unique. 
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